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Abstract 

The Fourier transform of the deeply virtual Compton scattering amplitude (DVCS) with respect 
to the skewness parameter ( = Q^/2p ■ q can be used to provide an image of the target hadron 
in the boost-invariant variable a, the coordinate conjugate to light-front time t = t + z/c. As an 
illustration, we construct a consistent covariant model of the DVCS amplitude and its associated 
generalized parton distributions using the quantum fluctuations of a fermion state at one loop in 
QED, thus providing a representation of the light-front wavefunctions of a lepton in a space. A 
consistent model for hadronic amplitudes can then be obtained by differentiating the light-front 
wavefunctions with respect to the bound-state mass. The resulting DVCS helicity amplitudes are 
evaluated as a function of a and the impact parameter b±, thus providing a light-front image of 
the target hadron in a frame-independent three-dimensional light-front coordinate space. Models 
for the LFWFs of hadrons in (3-1-1) dimensions displaying confinement at large distances and 
conformal symmetry at short distances have been obtained using the AdS/CFT method. We also 
compute the LFWFs in this model in invariant three dimensional coordinate space. We find that 
in the models studied, the Fourier transform of the DVCS amphtudes exhibit diffraction patterns. 
The results are analogous to the diffractive scattering of a wave in optics where the distribution in 
a measures the physical size of the scattering center in a one-dimensional system. 
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I. INTRODUCTION 



Deeply virtual Compton scattering (DVCS), 7*(g) + p{P) — ^ 7('?') + piP'), where 
the virtuality of the initial photon = —q"^ is large, provides a valuable probe of the 
elementary quark structure of the target proton near the light-cone. At leading twist, QCD 
factorization applies flj , and each DVCS helicity amplitude factorizes as a convolution in x 
of the hard 7*g Compton amplitude with a hadronic sub-amplitude constructed from 

the Generalized Parton Distributions (GPDs) H{x,(,t), E{x,(,t), H{x,(,t) and E{x,(,t)- 
Here x is the light cone momentum fraction of the struck quark; the skewness C = fp:^ 
measures the longitudinal momentum transfer in the DVCS process. 

Measurements of the momentum and spin dependence of the DVCS process in p — > 
7 p can provide a remarkable window to the QCD structure of hadrons at the amplitude 
level. The interference of the DVCS amplitude and the coherent Bethe-Heitler amplitude 
leads to an asymmetry which is related to the real part of the DVCS amplitude {2!. 
The imaginary part can also be accessed through various spin asymmetries p]. In the 
forward limit of zero momentum transfer, the GPDs reduce to ordinary parton distributions; 
on the other hand, the integration of GPDs over x reduces them to electromagnetic and 
gravitational form factors. 

The DVCS helicity amplitudes can be constructed in light-cone gauge from the overlap 
of the target hadron's light-front wavefunctions. . Since the DVCS process involves off- 
forward hadronic matrix elements of light-front bilocal currents, the overlaps are in general 
non-diagonal in particle number, unlike ordinary parton distributions. Thus in the case of 
GPDs, one requires not only the diagonal parton number conserving n n overlap of the 
initial and final light-front wavefunctions, but also an off-diagonal n + 1 ^ n — 1 overlap, 
where the parton number is decreased by two. Thus the GPDs measure hadron structure 
at the amplitude level in contrast to the probabilistic properties of parton distribution 
functions. 

The GPDs have become objects of much theoretical as well as experimental attention 
since they provide a rich source of information of hadron structure. Burkardt has noted 
that a Fourier Transform (FT) of the GPDs with respect to the transverse momentum 
transfer Aj^ in the idealized limit C = measures the impact parameter dependent parton 
distributions h±) defined from the absolute squares of the hadron's light-front wave 
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a a. 



functions (LFWFs) in x and impact space [g, [Zj. The impact representation on the hght- 
front was first introduced by Soper 8| in the context of the FT of the elastic form factor 
(see Appendix . The function q{x, bj_) is defined for a hadron with sharp plus momentum 
P+, localized in the transverse plane, such that the transverse center of momentum vanishes 
R± = 0. [ One can also work with a wave-packet localized in the transverse position space 
in order to avoid a state normalized to a 5-function. ] Thus q{x, b±) gives simultaneous 
information on the distributions of a quark as a function of the longitudinal light-front 
momentum fraction x = k^/P^ = {k^ + k^)/{P^ + P^) and the transverse distance b± of 
the parton from the center of the proton in the transverse plane. We use the standard LF 
coordinates P^ = P° ± P^, = ±y^. Since the proton is on-shell, P^P^ — P^ = M^. 

Since the incoming photon is space-like (g^ < 0) and the final photon is on-shell (g'^ = 0), 
the skewness ( is never zero in a physical experiment. In this paper, we will investigate 
the DVCS amplitude in the longitudinal position space by taking the FT with respect to 
(. We show that the FT of the DVCS amplitude in ( reveals the structure of a hadron 
target in a longitudinal impact parameter space. Thus, our work is suited for the direct 
analysis of experimental data and is complementary to the work of Burkardt and Soper. 
Physically, the FT of the DVCS amplitude allows one to measure the correlation within 
the hadron between the incoming and outgoing quark currents at transverse separation b± 
and longitudinal separation a = b^P^/2 at fixed light-front time t = z + t/c. Since Lorentz 
boosts are kinematical in the front form, the correlation determined in the three-dimensional 
6_L, (T space is frame-independent. ^ 

Even though light-front dynamics was proposed by Dirac [9] more than fifty years ago, 
and the utility of the light-front momentum fraction x = k~^/P~^ dates to the inception of 
the Feynman parton model, very little is known about the longitudinal coordinate space 
structure of hadron wave functions and related physical observables. To the best of our 
knowledge, the first work to investigate this subject is Ref. [l^ where it is shown that by 
Fourier transforming the form factors one observes profiles in b~ with kinks and anti-kinks. In 
addition to the light-front longitudinal structure of DVCS amplitudes in one-loop QED and 
meson models, we also present the corresponding structure of the light-front wavefunctions 
(LFWFs) of the quantum fiuctuations of a lepton to order in QED. 

Burkardt |6| has noted the possibility of taking the FT with respect to the longitudinal 
momentum of the active quark. However, since the GPDs depend on a sharp x, the Heisen- 
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berg uncertainty relation severely restricts the longitudinal position space interpretation of 
GPDs. In contrast, we will deal directly with DVCS amplitudes which are integrated over 
X and take the FT with respect to the longitudinal momentum transfer. 



It has been shown in llj that one can define a quantum mechanical Wigner distribution 
for the relativistic quarks and gluons inside the proton. Integrating over k~ and /c"*", one 
obtains a four dimensional quantum distribution which is a function of r and k'^ where r is 
the quark phase space position defined in the rest frame of the proton. These distributions 
are related to the FT of GPDs in the same frame. However, the Wigner distributions cannot 
be measured experimentally. 

In contrast, we will study the observable DVCS amphtudes directly in longitudinal po- 
sition space. We shall show that the Fourier transforms of the DVCS amplitudes in the 
variable cr = where the three-dimensional coordinate b = is conjugate to 

the momentum transfer A, provides a light-front image of the target hadron in a frame- 
independent three-dimensional light-front coordinate space. We find that in the models 
studied, the Fourier transform of the DVCS amplitudes exhibit diffraction patterns. The 
results are analogous to the diffractive scattering of a wave in optics where the distribution 
in a measures the physical size of the scattering center in a one-dimensional system. 

A summary of our main results has been given in 12|. In this paper we will present a 
detailed analysis and provide several additional results. 

In order to illustrate the general framework, we will present an explicit calculation of the 
FT of the DVCS on a fermion in QED at one-loop order jl3|. In effect, we shall represent a 
spin-| system as a composite of a spin-| fermion and a spin-1 vector boson, with arbitrary 
masses This one-loop model is self consistent since it has the correct interrelation of 
different Fock components of the state as given by the light-front eigenvalue equation fl^ . 
In particular, its two- and three-body Fock components can be obtained analytically from 
QED. This model has been used to calculate the spin and orbital angular momentum of 
a composite relativistic system 1^ as well as the GPDs in the impact parameter space 
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17(1 . The calculation is thus exact to 0{a), and it gives the Schwinger anomalous 



magnetic moment, the corresponding electron's Dirac and Pauli form factors jisl . [l6 | as 
well as the correct gravitational form factors, including the vanishing of the anomalous 
gravitomagnetic moment B{0) in agreement with the equivalence theorem 1^. In addition, 
it provides a template for the wave functions of an effective quark-diquark model of the 
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valence Fock state of the proton light-front wave function. 

Deep inelastic scattering structure functions and their connection to the spin and orbital 
angular momentum of the nucleon have been addressed for a dressed quark state in light- 



front QCD in 
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using a similar Fock space expansion of the state. This approach 
has also been used to investigate the twist-three GPDs in 2l|. We will also present here 
numerical results for a simulated model for a meson-like hadron, which we obtain by taking 
a derivative of the dressed electron LFWF with respect to the bound state mass M, thus 
improving the behavior of the wave function towards the end points in x. In this model, 
the DVCS amplitude is purely real. A similar power-law LFWF has been used in 22| to 
construct the GPDs for a meson. 

In principle, the LFWFs of hadrons in QCD can be computed using a nonperturbative 
method such as Discretized Light Cone Quantization (DLCQ) where the LF Hamiltonian 
is diagonalized on a free Fock basis This has been accomplished for simple confining 
quantum field theories such as QCD{1 + 1) j^ ]. 

Models for the LFWFs of hadrons in (3 + 1) dimensions displaying confinement at 
large distances and conformal symmetry at short distances have been obtained using the 
AdS/CFT method [24]. We will also present the LFWFs in this hadron model in invariant 
three dimensional coordinate space by Fourier transforming in both x and k-^. 

The plan of the paper is as follows: Section II summarizes the kinematics of the DVCS 
process. In Section III we give the analytic expressions for the DVCS amplitude and its 
explicit formulae in QED at one loop. The calculation of the Fourier transform of the DVCS 
amplitude for an electron target at one loop is given in Section IV. The simulated hadron 
model is discussed in Section V. We then derive the DVCS amplitudes using a model meson 
LFWF as obtained from holographic QCD in Section VI. A summary and the conclusions 
are given in Section VII. 



II. KINEMATICS 

The kinematics of the DVCS process has been given in detail in 0, 3]. One can work 
in a frame where the momenta of the initial and final proton has a A ^ —A symmetry 
^; however, in this frame, the kinematics in terms of the parton momenta becomes more 
complicated. Here, we shall use the frame of Ref. The momenta of the initial and final 
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proton are given by: 

p^(p^o-..|^;). (2.1) 

where M is the proton mass. The four-momentum transfer from the target is 

A = P-P' = (^Ci^+, A^, ^±^j , (2.3) 
where t — /S?'. In addition, overall energy-momentum conservation requires A~" = P" — P'~, 

— * ~ 

which connects Aj_, (, and t according to 

i = 2P-A = _CA4_+A^ _ ^2.4) 

1-C 

The coordinate b conjugate to A is defined by 6 • A = |&'''A~ -|- |&~A+ — b± • A±. We also 
define the boost invariant variable a — 6~P+/2 so that |fe~A+ = |6~P+C = a(. Thus cr is 
an 'impact parameter' but in the boost-invariant longitudinal coordinate space. 

It is convenient to choose a frame where the incident space-hke photon carries q'^ — so 
that q"^ = —Q^ = —q]^ (however, it is not mandatory to choose this frame): 

0,.x,^p^ + ^j3^j , (2.5) 

(cP-.fx + A,, . (2.6) 

We will be interested in deeply virtual Compton scattering, where is large compared to 
the masses and —t. Then, we have 

C (2.7) 



2P-q 

up to corrections in Thus C, plays the role of the Bjorken variable in deeply virtual 

Compton scattering. For a fixed value of —t, the allowed range of C, is given by 
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III. DEEPLY VIRTUAL COMPTON SCATTERING 



— > 

The virtual Compton amplitude Mf^'^{q±, A±, Q, i.e., the transition matrix element of 
the process 7*(g) + p{P) —>■ 7(5'') +p{P'), can be defined from the light-cone time-ordered 
product of currents 

M^^(g-1,A^,C) = t J d'ye-'^-y{P'\TJ^{y)r{0)\P) , (3.1) 

where the Lorentz indices fi and u denote the polarizations of the initial and final photons 
respectively. In the limit Q"^ ^ (yo at fixed ( and t the Compton amplitude is thus given by 



M^-^(g-l,A^,C) = 6;e:^M^^(g-l,A^,C) = '4 ^ 
X I t'J{zX)U{P') 



^ UiP) } , (3.2) 



where P = ^{P' + P). For simplicity we only consider one quark with fiavor q and electric 
charge Cg. We here consider the contribution of only the spin-independent GPDs H and 
E. Throughout our analysis we will assume the Born approximation to the photon-quark 
amplitude; i.e., the "handbag" approximation, corresponding to setting the Wilson line to 
1 in light-cone gauge. In principle there can be rescattering corrections in the light-cone 
gauge between the spectators at leading twist analogous to those which occur in diffractive 



deep inelastic scattering [25|, but these will not be considered here. 



For circularly polarized initial and final photons (J, J are f or J,)) we have 



z — ie z — ( + ie 
t^^{z,C) = t^\z,C) = . (3.3) 

The two photon polarization vectors in light-cone gauge are given by 

where k denotes the appropriate photon momentum. The polarization vectors satisfy the 
Lorentz condition k ■ e = 0. For a longitudinally polarized initial photon, the Compton 
amplitude is of order 1/Q and thus vanishes in the limit Q"^ 00. At order 1/Q there are 
several corrections to the simple structure in Eq. (13. 2p . We do not consider them here. 
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The generalized parton distributions H, E are defined through matrix elements of the 
bilinear vector and axial vector currents on the light-cone: 

dy- 



Stt 



zP+y-/2 
1 



y+=0,y_i_=0 



2P^ 



U{P') 



2M 



U{P) 



(3.5) 



The off- forward matrix elements given by Eq. fl3.5l) can be expressed in terms of overlaps of 

nn 

LFWFs of the state We now calculate the matrix elements in terms of the LFWFs. 

For this, we take the state to be an electron in QED at one loop and consider the LFWFs 
for this system. 



A. DVCS in QED at one Loop 



The light-front Fock state wavefunctions corresponding to the quantum fluctuations of a 
physical electron can be systematically evaluated in QED perturbation theory. The light- 
cone time ordered contribution for the state to the DVCS amplitude are given in Fig. 6 of 
[4]. The state is expanded in Fock space, giving contributions from | e~7) and | e~e~e"^), in 
addition to renormalizing the one-electron state. The two-particle state is expanded as. 



T 

two particle 







dx d^k 



+ - + 1 ; xP^ 



'x(l — x) IQn^ 



1 
'2 
1 
2 



1; xP^ 
1; xP^ 



k, 



(3.6) 



Here and denote 



■- +1; a;P+, k^^ + ^\_^{x,k^^ 

where the two-particle states |s£, s^; x, k_i_) are normalized as in P| 
the 2;-component of the spins of the constituent fermion and boson, respectively, and the 
variables x and k± refer to the momentum of the fermion. The light cone momentum fraction 
satisfy < Xj < 1, Yli^i — 1- employ the light-cone gauge = 0, so that 



Xi 



p+ 



the gauge boson polarizations are physical. The three-particle state has a similar expansion. 
Both the two- and three-particle Fock state components are given in [^]. The two-particle 
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wave function for spin-up electron are [4, Il3|, [1 



2 

V^Ii ,i(x,fcx) = -v^(M-f)^, 

2 

^Ii ,(x,fcx) = 0, 



(3.7) 



9?(a;, /cj 



X 



M2 



1-x 



(3.8) 



Similarly, the wave function 



Following references 



'or an electron with negative helicity can also be obtained. 



Il5l |. we work in a generalized form of QED by assigning a 
mass M to the external electrons, a distinct mass m for the internal electron lines, and a 
nonzero mass A for the internal photon lines, assuming the stability condition M < m + A. 
This provides a model for a composite fermion state with mass M with fermion and vector 
"diquark" constituents. The electron in QED also has a one-particle component 

dx d'^k 



one particle 







where the single- constituent wavefunction is given by 



1671^6(1 -x) 6\k^) 



±-:xP+,ki 



(3.9) 



(3.10) 



Here V Z is the wavefunction renormalization of the one-particle state and ensures overall 
probability conservation. Since we are working to 0{a), we can set Z = 1 in the 3-^1 
wavefunction overlap contributions. At x = 1, there are contributions from the overlap of 
one particle states which depend on Z. We have imposed a cutoff on x near this point. Also, 
in order to regulate the ultraviolet divergences, one has to introduce a regulator. Here, we 
use a cutoff A on the transverse momentum fc-*- as a regulator. 

In the domain ( < z < 1, there are diagonal 2^2 overlap contributions to Eq. 
fl33|) . both helicity flip, Ff_ (A' ^ A) and hehcity non-flip, Ff+ (A' = A) [4]. The GPDs 
H{2^2){z, (, t) and £'(2-^2) (-2^, C; t) are zero in the domain ^ — 1 < 2 < 0, which corresponds to 
emission and reabsorption of an from a physical electron. Contributions to H(^n->n){z, t) 
and E(^n^n){z,(,t) in that domain only appear beyond one-loop level since the DVCS am- 
plitude contains integrations over z, y~, and x. When the integration over y~ is performed, 
the fermion part of the bifocal current yields a factor 6{z — x) and the anti- fermion part of 
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the bilocal current yields a factor 5{z + x). The latter contribution is absent in the one- loop 
DVCS amplitude of a electron target, which we consider in the present work. 
We have, 



p22 



4(i-i)^w 



167r3 



^22 



1 {A^-iA' 
2M 

d2^ 



-E(2^2)(a;,C,^) 



where 



Using the explicit form of the two-particle wave functions, we obtain, 
"2 1 [{l-C)+x{x-0] 



p22 _ 



167r3 (1 - x) 



We use the notation 



(3.11) 



(3.13) 



/f ^ + 1^ + [B{x, C) + M'x{l -x)-m\l- x) 



M-—\ -^)(^-C) .ATF 



(3.14) 



L, = (fc^)^-2^^.A^|i-^-S(x,C), 
L2 = {k^f - M^x{l -x)+ rn^il - x) + }?x 



(3.15) 



and S(a;, C) = — ^jz^ 



2 (l-a:) _ y2 (a:-C) 



(1-C) 



rTVF 



rATF 



L1L2 



Trlog 



(1-0 



. The integrals are given by 



TT / df3 



—M'^x{l — x) + m?{l — x) + X^x 
1 



"i?(x,C,/3) 



(3.16) 



where L>(x, C, /3) = l3m'^{l-x)-l3M'^x{l-x)+l3X^x-{l-l3)B{x,C,)-{l-l3f{l-x'f{A^y^ 
Here A is the cutoff on the transverse momentum and x' = jfE^- 
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The helicity-flip part can be written as, 



n22 



x[x — 



m 



M 



+ 



m \ x{x — () 



7r3Vl-C C-xJ 



-'1 ) 



(3.17) 



where l[ = J ^^1^^ and I2 — I ~T^^^- have used the notation Ay = — iA^. These 
integrals can be done using the method described in Section III of Ref. jl?!] and we obtain 



=71 I dy 



y (l-x')A 



V 



(3.18) 



Q(a;,C,y) 

where Q{x,C,y) = (1 - y)[-M'^x{l - x) +m^{l - x) + X^x] - yB{x,() - 1/^(1 - x')^{A±f 



and 



(3.19) 



/o Qi.xX^y)' 

The scale A dependence is suppressed in F'^_. 

The contribution in the domain, < 2; < C comes from an overlap of three-particle and 
one-particle LFWFs. When the electron's helicity is not flipped, this contribution is given 
by B, 



n31 



^_ / i^(3^i)(x,C,t) - 



4(l-§)v/W 



^(3^1)(x,C,t) 



'1-C 



16^ 



^1 1 1 (a;, 1 - C,C - a:, -A_L, Aj 



"2 +3 "3 



+ , ^(x,l-C,C-a;,/i;±,-Ax,A_ 



"2 +7 

and can be written as, using the three-particle wave function, 
^2 



(3.20) 



C (1-a;) L(l-x) 



j/NF 
^1 



(1-0 



J/NF 
■^2 



e2 v^r^ x^ (C - x) 



M - 



m 



M 



m 



xJ\il-0 iC-x) 



T'NF. 



(3.21) 



with 



/NF 



iNF 



J: 



pNF 
"^3 



J^' + [M'x{l 
2x 



X) 



x] 



x^x].4 



jNF 
•^2 



jNF ^ ^M^^{1 - x) - m\l - x) - X^x + A{x, C)] J3 



jNF 
'^3 ■ 



(3.22) 
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We denote 



= {k^f - M^x{l -x)+ m^{l -x) + X^x, 



and C) 
written as: 



C(i-0 



h = {k^r - 2A;^ ■ + Mx, 0- (3.23) 
x{l - x)(A-L)2 + m2((l -() + x C (C - x)M2]. The integrals can be 



NF 



J. 



NF 



J. 



NF 



vrlog 
vrlog 



-M^x(l — x) + ni?{l — x) + X^x I- 

A2 



1 



(3.24) 



^2 io C'(x,C,/3)' 

where C(x, C, /3) = -pM^x{l - x) + pm^{l - x) + (3X^x + (1 - (3)A{x, C) - (1 - Pf^Al 
The hehcity-flip part is given by 
1 (Ai-iA2^ 



n31 



'1-C 



2M 
d^fc, 



E(3_l)(x,C,t) 



+7 +7 "7 



(3.25) 



Using the three-particle wave function, this can be written as, 



8n'^ \ X / 



m\ x^(C — x) 



c 



1 



(1-0 



7^ 

•^2 



1 



(l-x) 



F 
^1 



+ 



87r3 



M 



+ 



m 



x2(C -x) 



where 



Jo 



L(i-C) (C-x)J C (l-x) 

d^k^^h 



^1 



(3.26) 



•i^^ (l-^)A^f 
n I dy- ^ 



(3.27) 



^1^2 " 7o Pix,C,y) 
with P(x, C, ?/) = (!- y)A(x, C) - (1 - 2/)^f^(A^)2 + y[-M^x{l - x) + m^{l - x) + X^x]. 

We calculate the DVCS amplitude given by Eq. (13. 2p using the off-forward matrix ele- 
ments calculated above. The real and imaginary parts are calculated separately using the 
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prescription 



dx -—-Fix, = P dx -Fix, C) - tnFiC, ()■ (3.28) 

x-C + ie Jq x-C 

Here P denotes the principal value defined as 

P f dx^Fix,C) = \im \ f -J—Fix,C)+[ -^Fix,0] (3.29) 



where 

Fix,0 = 4^(x,C,A^), for 0<x<C 
= F^^ix,C,A^), forC<a;<l 

with ij = ++ for helicity non-flip and ij = H — for helicity flip amplitudes. Since the off- 
forward matrix elements are continuous at x = -^(C?C) = Ffj'ix = CX^^±) = Ff^ix = 
(,(,A±). Note that, for an electron state, the contribution vanishes for x < and the 
principal value prescription cannot be used at x = 0. The off-forward matrix elements F^^ 
(which contribute in the kinematical region < x < Q vanish result there 

is no logarithmic divergence at this point for nonzero (. But, we need to be careful here as 
when we consider the Fourier transform in a space, ( can go to zero and divergences from 
small X can occur from which is finite and nonzero at x, C — ^ 0. 

The imaginary part of the amplitude when the electron helicity is not flipped is then 
given by 

Im[M++](C, Ax) = ne'Fllix = C,C,^±)- (3.30) 

A similar expression holds in the case when the electron helicity is flipped (Im[M^ ](C, A-*-)) 

in which are replaced by The helicity- flip DVCS amplitude is proportional to 

(Ai — iA2) as seen from Eqs. (13.171) and fl3.26p . In the numerical results for the helicity 
flip processes that we present here, for simplicity we have taken A2 = 0. The imaginary 
part receives contributions from x = (. The off- forward matrix elements are continuous at 
X = (, and in Eq. (I3.30p we have used this continuity. The other regions of x contribute to 
the real part. It is to be emphasized that we are using the handbag approximation of the 
DVCS amplitude. Contributions from the Wilson lines are in general not zero, and they can 
give rise to new phase structures as seen in single-spin asymmetries [26|. 



14 



The real part of the DVCS amphtude is given by, 



Re [M++] (C,Aj 



-e" I dx F: 



dx F^i^(x,C,A_ 



1 1 
- + ^ 



1 

X 



X 



c 



(3.31) 



Similar expression holds for the helicity flip DVCS amplitude. 

The energy dependence of the DVCS amplitude at the high energies s » Q'^,—t can 
be deduced up to logarithms using Regge analysis. In our QED model for DVCS, two spin 
1/2 propagators are exchanged in the t-channel. This implies Regge behavior s""^^^ where 
a(0) = ji + j2 — 1 = 0. Thus, up to logarithms, the DVCS amplitude has no dependence on 
the initial energy of the incident electron and photon. 

IV. CALCULATION OF THE FOURIER TRANSFORM 



In order to obtain the DVCS amplitude in b space, we take a Fourier transform in ( as. 



2tx 



(4.1) 



where o = ^P^b~ is the boost invariant longitudinal distance on the light-cone. The spatial 
properties of deep inelastic scattering obtained from a Fourier transform of structure func- 
tions from X = k^/P^ to space has been discussed by Hoyer 27|. In Fig. [T]we show the 
handbag diagram of the DVCS amplitude in coordinate space, which is similar to Fig. 10 
of the above reference. 

Both real and imaginary parts of the DVCS amplitude are obtained separately. The real 
part of the amplitude depends on the cutoffs. Since the off-forward matrix elements are 
continuous at x = (, the DVCS amplitude is independent of the cutoff ei. The cutoffs have 
to be chosen such that €2 — 61 > e, €2 + ^1 < 1 — e, in order to have the correct principal value 
integration. In our numerical analysis, we have taken e = ei = 62/2 = 0.001. As stated 
before, the cutoff at x = is imposed for the numerical calculation and has a small effect 
on the result. If, instead of imposing a cutoff on transverse momentum. A, we imposed a 



cutoff on the invariant mass^lSj, then the divergences at x = 1 would have been regulated 
by a non-zero photon mass [28,]. The DVCS amplitude at a; = 1 also receives a contribution 
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FIG. 1: The handbag diagram for DVCS amphtude viewed in coordinate space. The position of 
the struck quark differs by x~ in the two wave functions (whereas = x-^ = 0). 



from the single particle sector of the Fock space [J, |16|, Il7|, l2l|] , which we did not take into 
account. A detailed discussion about the cutoff scheme is provided in Appendix [Bl 

All Fourier transforms have been performed by numerically calculating the Fourier sine 
and cosine transforms and then calculating the resultant by squaring them, adding and 
taking the square root, thereby yielding the Fourier Spectrum. The amplitude is divided 

4 

by the normalization factor In Fig. [21 we have shown the two particle LFWFs of the 
electron as a function of x for different k-^. We have taken m = 0.5 MeV, M = 0.51 MeV 
and A = 0.02 MeV. The wave functions are similar for a slight change of parameter values, 
however for m < M, there will be a node in tp-i/2+i{x, k^) at m = xM, which is seen in 
Fig. 1 (c). The effect of the node is almost negligible for these parameter values. 

In Fig. [3] we have shown the Fourier Spectrum (FS) of the 2-particle LFWFs given by 
Eqs. (13.71) . for the same mass parameters as in Fig. [21 The Fourier transform (FT) has been 
taken with respect to x for fixed values of transverse momentum k^. The wave functions 
r]i/2+i, V-i/2+1 and ?7i/2_i are obtained as. 
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FIG. 2: (Color online) Two-particle LFWFs of the electron vs. a; for M = 0.51 MeV, m = 0.5 
MeV, A = 0.02 MeV and fixed values of \ \= k in units of MeV. In (b) and (c) we have divided 
the LFWFS by the factors {k^ + ik^) and {—k^ + ifc^) respectively. 



where Xk 



is the peak of the distribution - where all the constituents in the 



n-particle Fock state have equal rapidity. 

All helicity components of the wave function show peaks at a = 0, the height of the peak 
sharply increases as /c"*" decreases and decays away from a = 0. In Fig. HI we have shown the 
complete Fourier transforms of the two particle LFWF for k=0.2 corresponding to the FS 
presented in Fig. [3]^b) and (c) to illustrate the difference. Though the real and imaginary 
parts of the FT (i.e., the cosine and sine transforms, respectively) individually exhibit a 
diffraction pattern, in Fig. |3](a) they are just out of phase to produce any diffraction pattern 
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FIG. 3: (Color online) Fourier spectrum of the two-particle LFWFs of the electron vs. a for 
M = 0.51 MeV, m = 0.5 MeV, A = 0.02 MeV and fixed values of | A;^ |= A; in MeV. In (b) and (c) 
we have divided the LFWFS by the factors {k^ + i/c^) and (—k^ + ifc^) respectively. 



in the FS. It is well-known in the theory of the Fourier representation of signals 29|] that the 
amplitude and phase play different roles and in some cases many of the important features 
of a signal are preserved only if the phase is retained. 

The plots of the DVCS amplitude have been done by fixing —t and varying both ( and 
Aj^. In Fig. [5](a) we have shown the imaginary part of the helicity flip DVCS amplitude 
M+_ as a function of ( for different values of —t. lm[M^_] is zero as C 0, increases 
continually with (, and then falls down sharply at the end. It increases for higher t for the 
same (. Fig. Mjo) shows the helicity non-flip part of the corresponding amplitude Im[M++] 
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FIG. 4: (Color online) Fourier transform (FT) of the two particle LFWFs of the electron vs. a for 
M = 0.51 MeV, m = 0.5 MeV, A = 0.02 MeV and fixed values of | /c^ |= in MeV. Re and Im 
denote the real and imaginary parts of the FT, and FS denotes the Fourier spectrum presented in 
Figs. Etb) and (c). 



vs. (. Unlike Im[M_|__], Im[M_|.+] is non- vanishing aX ( = 0, it decreases for higher —t. 
The largest allowed value of ( is given by Eq. ( 12. Sp for fixed t. Fig. [6](a) shows the plot 
of the real part of the helicity flip DVCS amplitude for the same values of the parameters. 

Re[M^ ] is non-vanishing at C — 0. For small (, it is almost fiat for a fixed —t and then 

falls down at large (. Fig. Mjo) shows the plot of the real part of the corresponding helicity 
non- flip amplitude M+_|_ vs. (. It shows a different functional behavior as it increases both 
for small as well as large (. 

Fig. [7] (a) shows the FS of the imaginary part of the helicity flip amplitude vs. a 
for M = 0.51 MeV, m = 0.5 MeV and A = 0.02 MeV. The peak of the FS of Im[M+_] 
(i.e.,|/m[y4+_] I) increases with | t \. The increasing behavior of the helicity flip amplitude 
7*6^(5"^ = 1/2) — > 76^ (S"! = —1/2) at small | t \ reflects the fact that one needs to transfer 
one unit of orbital angular momentum AL^ = ±1 to the electron to conserve J^. Note that 
the initial and final photon are taken to have transverse polarization. A similar behavior 
is expected for Yp\Sp = 1/2) -> -fp^iS^ = -1/2) amplitude. The FS of Im[M+_] does 
not show a diffraction pattern in a. In Fig. [7] (b) and (c) we have shown the FS of the 
imaginary part of the helicity non-flip DVCS amplitude vs. a for the same parameter values. 
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FIG. 5: (Color online) Imaginary part of the DVCS amplitude for an electron vs. C for different 
values oft : (a) helicity-flip part, (b) helicity non-flip part. We have taken M = 0.51 MeV, m = 0.5 
MeV, A = 0.02 MeV. The parameter t is given in MeV. 
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FIG. 6: (Color online) Real part of the DVCS amplitude for an electron vs. C for different values 
of t : (a) helicity-flip part, (b) helicity non-flip part. We have taken M = 0.51 MeV, m = 0.5 MeV, 
A = 0.02 MeV. The parameter t is given in MeV. 

The helicity non-flip amplitude depends on the scale A. We have taken A = Q. Fig. [7] (b) 
shows the plot for Q = 10 MeV, 5(c) is for Q = 50 MeV. Im shows a diffraction 

pattern, the peak becomes narrower as | t | increases. In Fig. |8]we have shown the complete 
Fourier transform of the imaginary part of the DVCS amplitude corresponding to the FS 
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FIG. 7: (Color online) Fourier spectrum of the imaginary part of the DVCS amplitude for an 
electron vs. a for different values of t : (a) when the electron helicity is flipped, (b) and (c) when 
the helicity is not flipped. In (b) Q = 10 MeV, in (c) Q = 50 MeV. The mass parameters are 
M = 0.51 MeV, m = 0.5 MeV, A = 0.02 MeV. The parameter t is given in MeV. 

presented in Fig. W^a). It again shows that the real (cosine transform) and the imaginary 
(sine transform) of the Fourier transform individually show a diffraction pattern, but they 
are out of phase and thus the FS does not show the diffraction pattern. 

The number of minima in the diffraction pattern increases with | t \ for fixed Q or, in 
other words, the first minima move in with increase of | t | . For Q = 50 MeV, the behavior is 
the same, the number of minima are higher for higher t for the same a range. The number 
and the positions of the minima are independent of Q. Only the magnitude of the peak 
changes with Q. Some of the plots of the FS of the DVCS amplitude show similarities with 
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FIG. 8: (Color online) Fourier transform (FT) of the imaginary part of the helicity flip DVCS 
amplitude for an electron vs. a for t=-5.0. Re and Im denote the real and imaginary parts of the 
FT. The mass parameters are M = 0.51 MeV, m = 0.5 MeV, A=0.02 MeV. 

the FS of the LFWFs themselves. The generalized parton distributions are related to the 
form factors, and the form factors can be written as overlaps of LFWFs. In fact, for a meson 
in 1 + 1 dimensional QCD, the form factor becomes an overlap of the LFWFs with different 
longitudinal momentum fractions, x [s^, and the contribution is similar to the 2^2 part 
of DVCS amplitude. 

Fig. [9] (a) shows the FS of the real part of the helicity flip amplitude vs. a, where we 
chose the same values of the parameters, M = 0.51 MeV, m = 0.5 MeV and A = 0.02 MeV. 
The FS i.e., |i?e[yl_|__]| shows a diffraction pattern in a. With the increase of —t, the central 
peak increases and its width decreases . The helicity flip part of the DVCS amplitude does 
not depend on Q. 

The non-existence of the diffraction pattern in the FS of the imaginary part of the helicity 
flip amplitude in a is due to its different behavior in (, as seen from Figs. [5]and[6l Im[M^_] 
decreases smoothly as ( decreases to vanish a.t ( = which is distinctly different from all 
other amplitudes including i?e[M+_]. All other amplitudes show some flatness or a plateau 
in ( and their FS in a space shows a diffraction pattern. 

In Fig. [9] (b) and (c), we have plotted the FS of the real part of the helicity non-flip 
DVCS amplitude vs. a. One can see the diffraction pattern here as well. Fig. Mjo) is 
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FIG. 9: (Color online) Fourier spectrum of the real part of the DVCS amplitude for an electron vs. 
a for different values of t : (a) when the electron helicity is flipped, (b) and (c) when the helicity 
is not flipped. In (b) Q = 10 MeV, in (c) Q = 50 MeV. The mass parameters are M = 0.51 MeV, 
m = 0.5 MeV, A = 0.02 MeV. The parameter t is given in MeV. 

for Q = 10 MeV and (c) is for Q = 50 MeV. As before, the qualitative behaviors of the 
diffraction pattern do not change with Q. For the same | t |, the number of minima and 
their positions are independent of Q for any fixed | t |, only the height changes with Q. For 
each Q, the peak at cr = is sharper and higher as | t | increases. 

Instead of using the C, and t variables, we can define another set of variables C, and T, 
where T is defined as T = j ■ The arguments of the final state LFWF then are x' = 
and k'^ = fc-*- — (1 — x)^/T, in other words, the transverse momenta become decoupled from 
(. We can now take ( and T as independent variables and the GPDs as well as the DVCS 
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FIG. 10: (Color online) Imaginary part of the helicity non-flip DVCS amplitude for an electron vs. 
a for fixed T in MeV . 



amplitude can be expressed in terms of them. They are however, connected through 



t 



c 



:i - c)T. 



(4.3) 



This relation determines the range of allowed values of C, and T, such that t « . In 
practice, C, can never become very close to 1. The C, dependence of the DVCS amplitude now 
comes purely from x' for fixed T. Fig. [10] shows a plot of the imaginary part of the helicity 
non-flip amplitude for fixed T vs. a. For very small cr, the slope of the a distribution is 
given by. 



da 



l-e 



(x, X, T) + -FY+(a;, x, T)). 



(4.4) 



Thus the slope and therefore the width of the a distribution depends on the second moment 
of the GPDs at X = C- 

V. SIMULATED BOUND STATES 

For the dressed electron state, the real part of the DVCS amplitude depends on the 
cutoffs at X = and x = 1. We have chosen the cutoff scheme discussed earlier. The cutoff 
at X = is taken for the numerical calculation and its effect on the result is small. However, 
starting from this QED point-like model where the electron fluctuates to spin-half plus spin 
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FIG. 11: (Color online) Fourier spectrum of the LFWFs of the simulated hadron model vs. cr for 
M = 0.5 MeV, m = A = 1.0 MeV and fixed values of | k-^ \= k in MeV. We have divided the 
LFWFs by the normalization constant. In (b) and (c) we have divided by the factors {k^ + ik'^) 
and {—k^ + ifc^) respectively as well. 

one constituents, one can construct LFWFs for the hadrons. In the two- and three-particle 
LFWFs for the electron, the bound state mass M appears in the energy denominators. A 
differentiation of the QED LFWFs with respect to improves the convergence at the end 
points: x = 0, 1, as well as at high k"^, thus simulating a bound state valence wavefunction. 
Differentiating once with respect to will generate a meson-like behavior of the LFWF. 
Thus we write the hadron two-particle LFWFs as 

d 

ips^s^ix, k±) = M^——iJs,s2{x, k±) (5.1) 



25 



where V'sisal^; ^-l) are the electron LFWFs. Taking the Fourier transform in k± also we can 
write the wavefunctions in a and transverse impact parameter b± as 

Xs.sM b±) = / c?'A;xe-(--^'=)e^^-^-V;,,,,(x, k^), (5.2) 

where the peak of the distribution Xk — / o '^o ■ Writing k± ■ b± — kb cos 9, where 
b = \b±\ and k — \k±\ and performing the integration over 9 we obtain 

Xs,s,{a, b^) = 2' dx J kdke"^^^-'>'\27r)Mbk)i>s,s,{x, k^), (5.3) 

where Jo{bk) is the Bessel function. For the wavefunctions Xi/2±i we have the explicit 
momentum components {ki ± ik2) present in the numerators. We use 

ik±-b± 



J d^k^ik.^ik^y^--'- = {-t) I d'k^{^^±t^y 



= 27ri J edk^^^^^Mkb) (5.4) 
For the plot of the wavefunctions we set 62 = in the above expression, 

/ ,N e /7r,,2. /■ , /■,2„ x2(l -a;)i/Vi(A;6)e*'^(^-*'=) 

(5.5) 

For computational purpose, we use 

1 Z"' 

J„(A;6) = - / d9cos{n9 - kbsm9). (5.6) 
TT Jo 



This procedure does not provide an actual model for a 'meson' wavefunction since the 
two constituents have spin-half and spin-one. However, if we differentiate once more we can 
simulate the fall-off at short distances which matches the fall-off wavefunction of a baryon. 
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FIG. 12: (Color online) Fourier Spectrum of the LFWFs of the simulated hadron model plotted 
in a, \b±\ space for M = 150 MeV, m = X = 300 MeV. X1jX2 and xs are | Xi/2+1 I; I Xi/2-1 I and 
I X-i/2+1 I) respectively. In the plots |6_l| runs from 0.001 to 0.01 MeV^^ and a runs from -25 to 
+25. 

in the sense that the form factor Fi(Q^) computed from the Drell-Yan-West formula will 
fall-off like In this case, we obtain a quark plus spin-one diquark model of a baryon. 
Convolution of these wavefunctions in the same way as we have done for the dressed electron 
wavefunctions will simulate the corresponding DVCS amplitudes for bound-state hadrons. 
Note that the differentiation of the single-particle LFWF will give a vanishing result and 
as a result, the 3^1 contribution to the DVCS amplitude vanishes in this model. The 
resulting 7*^ ^yp DVCS amplitude has both real 2] and imaginary parts Q|. 

If we consider a dressed electron, the imaginary part from the pole at x = ( survives 



27 



0.005 



0.004 



0.003 



0.001 



0.0015 



(a) 



+ 




- V 


-t) 


=100 


Oi 




- V 


-t) 


=316 


0.002 




... V 


-t) 


=707 



0.2 0.4 0.6 0.8 



0.001 



\ I 
\ I 

\ 
\ 



0.0005 



V(-t)=ioo 

V(-t)=316 
V(-t)=707 




FIG. 13: (Color online) Real part of the DVCS amplitude for the simulated hadron state. The 
parameters are M = 150, m = X = 300 MeV. (a) Helicity flip amplitude vs. (b) Fourier spectrum 
of the same vs. a. The parameter t is in MeV^. 

because of the numerator factor in the electron's LFWF. This numerator behavior 
reflects the spin-1 nature of the constituent boson. The x — —>■ singularity is shielded 
when we differentiate the final state n = 2 and n = 3 LFWFs with respect to and, 
as a result, the imaginary part of the amplitude vanishes in this model. We thus have 
constructed a model where the DVCS amplitude is purely real. However, the forward virtual 
Compton amplitude 7*p '~f*p (whose imaginary part gives the structure function) does 
not have this property. The pole at x = ( is not shielded since the initial and final n = 
2 LFWFs are functions of x. It is worthwhile to point out that in general the LFWFs 
for a hadron may be non- vanishing at the end points [31], and recent measurements of 
single spin asymmetries suggest that the GPDs are non- vanishing at x = ( 32(]. A more 



realistic estimate would require non- valence Fock states [33|]. An equivalent but easier way 
to construct the hadronic model is to differentiate the DVCS amplitudes with respect to the 
invariant mass squared (M^) of the initial and final bound states. Thus one can calculate 
the quantity Mp-^^Mf-^jjpAij{Mj, Mp) where M/, Mp are the initial and final bound state 
masses. 



For numerical computation we use the discrete (in the sense that the denominator is 
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FIG. 14: (Color online) Real part of the DVCS amplitude for the simulated hadron bound state. 
The parameters are M = 150, m = A = 300 MeV. (a) Helicity non-flip amplitude vs. C, (b) Fourier 
spectrum of the same vs. a, (c) Structure function vs. x. The parameter t is in MeV^. 



finite) version of the differentiation: 



-, A{Mf)-A{Ml) 
6M^ 



(5.7) 



where = (Mf + M|)/2 and 6M^ = {Mf - M|). Then we have 



d 



d 



MfM^p 



(5.8) 



Aij{Mi, Mp) is the DVCS amphtude for an electron. The differentiation with respect to Mp 
of the amplitudes (helicity flip and non-flip) brings in an extra factor oi x — C, and thus the 
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imaginary parts of the DVCS amplitudes vanish in this model, as discussed before. The real 
parts of the DVCS amplitudes survive and show diffraction patterns. We take Mii,Mfi = 
150 + 1.0 and Mj2, Mf2 = 150 - 1.0 and the fixed parameters M = 150, m = A = 300 MeV. 
In Figs. [TT] and [T2] we have shown the FS of the 2-particle LFWFs for this model. The 
parameter values are scaled in Fig. [TT], but the qualitative behaviours are the same. Since 
the wave function now vanishes at x = 0, 1, the FS is localized, and it decays sharply beyond 
I a 1= 10. The peak decreases more sharply for higher k± or lower In Fig. [13] (a) and 
(b) we have shown the helicity-flip DVCS amplitude for the hadron model as functions of ( 
and its FS as a function of a respectively. In Fig. [14] (a) and (b) we have shown the real 
part of the helicity non-flip amplitudes of the same model as a function of ( and its FS as a 
function of a respectively. Notice that the helicity non-flip part of the amplitude no longer 
depends on the scale. The amplitude decreases as ( increases, in contrast to the behavior 
for the electron. The FS of both the helicity-fiip and non-flip DVCS amplitudes show a 
diffraction pattern in a. Fig. [T4](c) illustrates the structure function F2{x) for this model as 
a function of x. 



A. An Optical Analog 

We propose an optics analog of the behavior of the FT of DVCS amplitude in a. The 
similarity of optics and quantum fields on the light cone was first explored long ago by 



hey established the similarity of paraxial-wave optics 



Sudarshan, Simon and Mukunda. 
and light cone dynamics of scalar 3j] and Maxwell equations |35|. In our case, we are 
effectively looking at the interference between the initial and final waves of the scattered 
proton. The final-state proton wavefunction is modified relative to the incident proton 
wavefunction because of the momentum transferred to the quark in the hard Compton 
scattering. The change in quark momentum along the longitudinal direction can be Fourier 
transformed to a shift in the light-front position of the struck quark; thus one can simulate 
a change in the quark's longitudinal LF coordinate by an amount a = . This is 

analogous to diffractive scattering of a wave in optics where a plays the role of the physical 
size of the scattering center in a one-dimensional system. We are using t to register the 
change in the transverse momentum of the quark in the scattering. The positions of the 
first minima move in with increasing \t\. 



30 



Notice that the integrals over x and C, are of finite range. More importantly, the upper 
limit of C integral is C,max which in turn is determined by the value of — t. The finiteness of 
slit width is a necessary condition for the occurrence of diffraction pattern in optics. Thus 
when the integration is performed over the range to Cmaxi the finite range acts as a slit of 
finite width and provides a necessary condition for the occurrence of diffraction pattern in 
the Fourier transform of the DVCS amplitude. 

When a diffraction pattern is produced, in analogy with single slit diffraction, we expect 
the position of the first minimum to be inversely proportional to Cmax- Since C,max increases 
with — t, we expect the position of the first minimum to move to a smaller value of a, 
in analogy with optical diffraction. In the case of the Fourier Spectrum of DVCS on the 
quantum fluctuations of a lepton target in QED, and also in the corresponding hadronic 
model, one sees that the diffractive patterns in a sharpen and the positions of the flrst minima 
typically move in with increasing momentum transfer. Thus the invariant longitudinal size of 
the parton distribution becomes longer and the shape of the conjugate light-cone momentum 
distribution becomes narrower with increasing \t\. 
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TABLE I: Positions of minima in the diffraction pattern for different —t for the simulated hadron 
model DVCS amplitudes. 

From Table [H we can see that for a flxed (— t), higher minima appear at the positions 
integral multiple of lowest minimum. This is consistent with the single slit diffraction law 
for n^^ minima: sin^^„ = nX/w where A is the wavelength of the light and w is the slit width. 
Here a plays the role of sin^„ (for large separation(L) between the slit and the detector in 
a single slit experiment, one can write sin^„ ^ ^) and the ratio X/w is determined by the 
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value of —t. Positions of the minima do not depend on the hehcity. The minima appear at 
the same places for both helicity flip and non-flip processes. 

We also observe a relation between the invariant momentum transfer squared —t and the 
distribution in a: the first minimum in a diffraction pattern is determined by 



sin 6*1 



_ A 

■-2 W 



(5.9) 



where ai is the position of the first minimum measured from the center of the diffraction 
pattern. Introducing another parameter to, with —t > —to we write the right hand side i.e., 
the ratio ^ as, where we have chosen —to = 2 x 10^ MeV^. Once to is fixed, the 

other parameter L can be found; ai = 8.5 MeV~^ for — t = 5 x 10'^ MeV^ gives 



8.5 



1 



0.2234 



/PTW log(3 X 104) 
which gives ~ 1376 MeV^^. Using this value of L we can compare with the other data: 





0-1 




1 




log(-t+to) 


264 


8.0 


0.211 


0.213 


316 


7.5 


0.198 


0.204 


707 


7.0 


0.185 


0.176 



TABLE II: Comparison of our proposed formula with the data, 
are in MeV~^. 



t is in MeV and the lengths 



Table HTl shows that our parameterization of ^ in terms of —t is quite accurate. 
If one Fourier transforms the change in transverse momentum to impact space b± 
6|, lZI], then one would have the analog of a three-dimensional scattering center. In this sense, 
studying the FT of the DVCS is very much like studying the Lorentz-invariant optics of the 
proton. In our analysis we have computed DVCS on an electron at 0{a) and its FT. Thus 
we have obtained the diffractive optics of the quantum fluctuations of an electron. 

It is interesting to recall that the scattering amplitude corresponding to an absorptive 
(i.e., negative imaginary) potential which is conflned to a sphere of flnite radius exhibits a 
diffraction pattern. For a classic treatment, see R. J. Glauber's lectures 36|]. For another 
example of diffraction patterns in the angular distribution of elastic proton-nucleus scattering 



using a multiple-scattering approach, see Ref. 

32 



37l |. It is worthwhile to remember that in 



FIG. 15: (Color online) The ground state ( L = 0, = 1) of two parton holographic light front 
wave function in 3D space. We have taken Aqcd = 0.32 GeV. Here runs from 0.001 to 6.0 
GeV"^ and a from -25 to 25. 

these examples the optical potential is complex. In our case, we observe diffraction patterns 
when we perform Fourier Transforms of real functions. 



VI. LFWF FOR MESON IN HOLOGRAPHIC QCD AND THE DVCS AMPLI- 
TUDE 

The normalized holographic QCD LFWF for the meson (qq) from AdS/CFT derived in 
ref. Q is 

^L,fe(a;, 6±) = BL,k\/^4^^^M^pL,kAQCD) (6.1) 

where BL,k = ^QCD[i-^)^TTJi+Lif3L,k)Ji-Lil3L,k)]~^^'^, ^ = A/a;(l - and (3L,k is the 

k-th zero of Bessel function J^. For ground state L = 0, k = 1 and we have 

0(x, b,) = v^o,i(x, b,) = AQc.y^(T^ ^°^f°'f^^^"^ (6.2) 



The corresponding momentum space LFWF is 2J] 



^(x,/€x) = V47r2 / e-'^'^-^^ 0(a;,&i±) . (6.3) 



In Fig. [151 we have plotted the two-parton bound state holographic LFWF from the 
AdS/CFT correspondence in 3D coordinate space after taking the FT of Eq. (16.21) with 
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respect to x. The ADS/CFT correspondence gives only the wavefunction for the qq sector. 
So, when we consider the DVCS amphtude with this wavefunction we can have contribution 
only from the 2 — * 2 process. Then the DVCS amplitude can be written as 

-1 



M(Aj_,C) 



E 



dx 



1 



X — C + ie X — ie 
1 — X 



d K± 



-A 



i-C ^ i-C 

1 — X 1 — X 



(6.4) 



C 1-C 

The transverse Fourier Transform of the DVCS amplitude gives the DVCS amplitude in 
the transverse impact parameter space b±. 



(6.5) 



Then 



A{b^X) = (2vr)^ / dxF{x,0 
+ 



1 — X 
1-x 1-C 



b±) (p{x, b 



1-C'l-x 



1 — X 



1-C i-x 



1-C 

&±)0(1 - a;, 6j 

1 — X 



(6.6) 



where -F(x, C) = ~^q{^ + Taking the FT of this DVCS amplitude with respect to C we 
obtain the amplitude in the 3-dimensional impact parameter space a, 6^. Substituting the 
wavefunctions given in Eq. ( 16.20 . we obtain 



A{a,b^_) 



1 

2^ 



2(27r) 



dce<~A{bU) 

4 ^QCD f 



rfCe^^^y dx F(x,C)Vx(x-C)[Jo(Xi)Jo(X2)](6.7) 



where 



= ^a;(l-x)^^|6x| /5o,iA 
1 — X 

X2 = ^{l-x){x-0-^\b^\ f3o,lAQCD ■ 

1 — X 

In Fig. [in] we show the FS of the DVCS amplitude in cr space for different fixed values of 
Again we see the diffraction pattern in the a space. 
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FIG. 16: (Color online) (a) The DVCS amplitude vs. ( and (b) the Fourier spectrum of the DVCS 
amplitude in the a space using the light front wave function for meson obtained from holographic 
QCD We have taken Aqcd = 0.32 GeV. Plots are in unit of e^. 6"*" is given in units of GeV~^. 

VII. SUMMARY AND CONCLUSIONS 

The deeply virtual Compton scattering process 7*p 7p provides a direct window into 
hadron substructure which goes well beyond inclusive measurements. The DVCS amplitude 
factorizes into the convolution of a hard perturbative amplitude, corresponding to Compton 
scattering on a quark current, with the initial and final state light-front wavefunctions of 
the target hadron. The LFWFs provide a general frame-independent representation of 
relativistic composite hadrons, and they are universal and process independent. 

In this paper, we have shown that the Fourier transform of the DVCS amplitude with 
respect to the skewness variable ( gives information of the proton structure in longitudinal 
impact parameter space a = ^P^b^ . 

As an illustration of our the general framework, we have worked with a simple relativistic 
spin 1/2 system, namely the quantum fluctuations of a lepton at one loop in QED. The 
different Fock components of the LFWFs in this case can be obtained from perturbation 
theory. Our calculation is exact to 0{a). This one-loop model provides a transparent basis 
for understanding the structure of more general bound-state systems. By differentiating the 
wavefunctions for the electron with respect to the square of the bound state mass M^, we 
have simulated bound state valence wavefunctions. 
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We have noted that there are two different types of overlaps contributing to DVCS when 
( is nonzero, namely a parton number conserving diagonal 2 —>■ 2 overlap and a 3 ^ 1 
overlap when an electron-positron pair of the initial state is annihilated. In fact, both these 
contributions are necessary in order to obtain ( independent form factors by taking the x 
moment of the GPDs. This invariance is due to the Lorentz frame independence of the light- 
front Fock representations of space-like local operator matrix elements and it reflects the 
underlying connections of Fock states with different parton numbers implied by the equation 
of motion. The Fourier transform of the amplitude with respect to ( involves both type of 
contributions in different kinematical regions. 

We have introduced the light-cone longitudinal distance a = and have shown the 

a dependence of the real and imaginary part of the DVCS amplitude. The DVCS amplitude 
in cr space represents an interference of the initial and final state LFWFs. We have also 
shown the a dependence of the LFWFs themselves. 

We have exhibited the light-front coordinate space structure of our model wavefunctions 
by performing the Fourier transform in the longitudinal and transverse momentum space. 
The wavefunctions exhibit diffraction patterns in the longitudinal coordinate space. We have 
presented the FS of the real part of the DVCS amplitude as well as the structure function 
in the models. The corresponding imaginary part of the model DVCS amplitudes vanish. 

Very recently, valence parton bound state holographic LFWFs from the AdS/CFT cor- 
respondence have been given 2J]. We have presented these wavefunctions for the meson in 
full three dimensional light-front coordinate space. We have also calculated the real part of 
the DVCS amplitude in the holographic model in light front longitudinal space for specific 
choices of the impact parameter {bj_). Again one observes diffraction patterns. Note that 
the imaginary part of the DVCS amplitudes vanish also in this model. 

Our analysis is the first to examine the longitudinal light cone coordinate a = h~P^ /2 
dependence of LFWFs and DVCS amplitudes. Our results for the DVCS amplitude in a are 
analogous to diffractive scattering of a wave in optics where the a distribution senses the 
size of the one-dimensional scattering center. Thus studying DVCS 7*p 7p in light-front 
longitudinal coordinate space is very much like studying the Lorentz-invariant optics of the 
proton. 
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APPENDIX A: RELATION BETWEEN BURKARDT AND SOPER DENSITIES 

The off-forward parton distribution appropriate for a spin-zero meson in the valence 
approximation is defined by 

' dz 5{x - z)e{z -of d^'- ~ ^ ^ ~ ^ 





dz 6{x - z)e{z - c) y d'h^ri^z^, - Y^^±)^i^^ ■ (Al) 



For skewness C = 0? 

H{x, C = 0, t) = H{x, A^) = j d^k^tp*{x, k±-{l- x)A^)i;{x, k^) 



(27r)2 

where the FT of the wavefunction is defined as 

i^{x,k^) = J ^e-^'--'-<P{x,b^) . (A3) 

The transverse Fourier transform of the zero skewness off-forward parton distribution H 
yields the impact parameter density function 

= ^^Y^, (A4) 

n 

where p{x,b±) is the Soper density defined in Eq. (5) of Ref. [8|. We thus find that the 
density obtained by Burkardt is the same as the Soper density. 
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APPENDIX B: REGULATORS 



Let us consider the real part of the DVCS amphtude for electron in one-loop QED given 



by 



Jo 



1 1 

- + 7 

X X — Q 

1 1 

- + 7 

X X — Q 



(Bl) 



which results after performing the transverse momentum integration. As described in the 
text, we use an ultraviolet cutoff A on the transverse momentum. 

The integrands for various DVCS amplitudes may exhibit singular behavior when x is 
near the end points. We also have a potential singularity when x ^ C, which can be regulated 
by using the principal value prescription. In the numerical work we implement the principal 
value prescription by employing suitable regulators and ensuring regulator independence in 
the limit where the regulator vanishes. We also note that we eventually integrate over C, 
which ranges between Cmm, which is close to zero and Cmax, which is determined by —t and 
approaches 1 in the limit —t^oo. 

The light-cone momentum fractions must remain positive. Let Cmm = C + for the 
second integral. Since Cmax = 1 — e, to make sure that x remains greater than C, in the second 
integral, we choose Xmax in the second integral to be 1 — e/2. Thus the regulated integral is 

K-e/2 



Re M{C,A±) 



e/2 



dx F-''{xX,A±) 



dx F'\x,C,A±) 



1 


+ 


1 




x 


X — 


c 


1 


+ 


1 




x 


X — 


c 



(B2) 



It is important to note that when ( is small, significant contribution to the integral comes 
from the second term that involves and when ( is large, significant contribution to the 
integral comes from F^^. For the helicity-flip case, both and F^^ vanish a.s x —* 0. For 
the helicity non-flip amplitude, F^^ vanishes as x ^ 0, but is finite clS X ^ 0. Thus 
the only potential problem at small x occurs in Eq. (lB2p for the second term involving 
when ( is small and we obtain a logarithmic divergence due to this problem. This is directly 
related to the non-vanishing of the electron wave function as x 0. For helicity non-flip, 
the function F also diverg 1. In this region, the DVCS amplitude also receives 

contributions from the single-particle sector of the Fock space which we do not take into 
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account in the present calculation. If one uses an invariant mass cutoff, the divergences 
at X — and x — 1 would have been regulated by non-zero electron and photon masses 
respectively. These regulators are not mandatory in our present calculations, and we have 
employed the simpler regulators as described above. 

APPENDIX C: AN ILLUSTRATIVE MODEL 

Let us start from the expression for the Fourier transform 

A{a) = ^JdC e^'^MiC) . (CI) 

The function A{a) is in coordinate space and the function M[C,) is in momentum space. 
Let us approximate M by the following function: 



M(C) = Mo for < C < Cma. 

= for C > Cma. (C2) 

One can use this step function to approximate a DVCS amplitude in which the dependence 
of the DVCS amplitude in C, is almost flat; in such a case we can take Mq = [M{C, — 

0)+M(C = Cma.]/2. 

We have. 



Mn 



max 



^ MoCmax sin(aCmax-/2) ^i^^^„^/2_ ^^^^ 

In this case, note that the cosine and sine transforms are completely in phase and the phase 
of the Fourier transform does not contain any extra information. The amplitude (i.e., the 
Fourier Spectrum) is given by 

I X I MoCmax I sin((TCmax/2) I , . 

I = aC»^/2 ■ 

The magnitude of the peak of the diffraction pattern 

A[a)max = — 1^ — (C5) 
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and the first diffraction minimum occurs at 

ai = . (C6) 

Smax 

In the case of the DVCS amplitude whose functional dependence on ( is very weak, we 
can further predict the position of the minima as follows. The extension of the function (max 
is given by 

Thus we find a precise relation between the minima of the diffraction pattern and —t. Since 
0"! is inversely proportional to (max which in turn increases with —t, the inward movement of 
the first minimum with increasing —t is readily explained. On the other hand the peak height 
is a product of Mq and (max- The amplitude Mq decreases monotonically with increasing 
— t. On the other hand (max increases with increasing —t. Thus the peak height has non- 
monotonic behavior with respect to —t. 

In Table Hm we compare the numbers with the case of the real part of the helicity non-flip 
amplitude presented in Figs. [TWa) and (b). Similarly one can also obtain estimates for the 
helicity flip amplitudes. 





Cmax 


Mo 


Peak (^^oCma.) 


First Minimum (ai = ) 


100 


0.48 


0.04 


3 X 10"^ 


13.09 


316 


0.84 


0.03 


4 X 10-^ 


7.48 


707 


0.96 


0.0175 


2.7 X 10^3 


6.54 



TABLE III: Simplified approximation for the real part of the helicity non-flip DVCS amplitude 
shown in Figs. [U] (a) and (b). The energy quantities are given in MeV and lengths in MeV~^. 

The essential ingredients for the diffraction pattern in the Fourier Spectrum are two 
characteristics of the DVCS amplitudes in the variable (: 

• A step (i.e., a sharp rise) and 

• a plateau. 

These are the essential characteristics of a function which is almost a constant that seems 
to be shared by the DVCS amplitudes that produce a diffraction pattern in the FS. The 
imaginary part of the helicity-flip DVCS amplitude for the electron state (Fig. [5]^a)) lacks 
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these properties and we do not observe any diffraction pattern in the corresponding Fourier 
Spectrum (Fig. [7|(a)). 

It is interesting to note that the simple model we have discussed in this appendix appears 
in antenna theory [38]. In the case of an aperture for which a uniform electric field is 
maintained over a finite distance, outside of which the field is zero, the angular spectrum 
which is the Fourier Spectrum of the aperture field distribution exhibits the diffraction 
pattern discussed in this appendix. 



APPENDIX D: DVCS AMPLITUDE IN THREE DIMENSIONS 

The significance of the amplitude in the boost-invariant a space can also be explained 
in the following way. The Dirac and Pauli form factors Fi{t) and F2{t) respectively, can be 
expressed in terms of the helicity non-flip part of the off- forward matrix element, 

^et us consider the dressed electron in the frame C = 0. The form factor can be written 

as y 

F{t) = C dx f d%^e-'^^-^^\ Mx,b^) f (Dl) 



The LFWFs in the mixed representation x,b^ are given by Eq. (1) of [24]. Note that the 
LFWFs are zero outside the region < x < 1. We denote : 

^2(x, b^) = iP2{x, b^)e{x)e{l - x) 

^3(2:1, X2, bi, b^) = i^sixi, X2, bi, bi)e{xi)e{x2)e{i - xi)e{i - X2) (D2) 

for the two particle LFWF and similarly for the three-particle wave function. We take 
FT of the LFWFs with respect to x ; and define 

/+00 
daie-''^^^'}^4x^,bi) (D3) 
-00 

where cxj are the boost invariant longitudinal distance on the light cone, conjugate to Xi = 
kf / . There are n — 1 independent cTj as well as b^. In terms of these, we can write, 

F{t) = J dx J dVe-'^^ '^y daij da2e*"^"e-'"^"$;(ai, 6^)<l>2((72, 6^) 

= 27r /" d%^e-'^^-'''^ [ da\ $2(a,6^) f (D4) 
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Note that as $„ are the FT of the wave functions \l/„ rather than ipn, it is mathematically 
correct to take the x- integrals from — oo to +00. When ( is non-zero, the form factor 
receives contributions from 2 — 2 and 3 — 1 components of the GPDs H and E. They can 
be obtained from, 



dx [F^iix, c, mx - c) + F^'4x, c, mc - x)] 



^/i^^plix, 1-cx-x, -b^)eic - x) 



1 — X 1 — X (1 — xy 



dx I d^^e-'^^-'"^ 



(D5) 



dx [FiHx, c, mx - c) + Fi'_{x, c, mc - x)] 

1 (Ai-^A2)(l-C/2) 



dx I d%^e-'^"-'" 



2M 



+r2\x\ --^)i^i{x, --^)^—^e{x - c) 
1 — x' 1 — x' (1 — xy 



1-0^ 



X' 



1-C" 



Here we have suppressed the explicit helicity indices and used 



' H{xX,t) 

' E{x,C,t) 
^ 2 



(D6) 



(D7) 



The form factors Fi{t) and F2{t) can be obtained in terms of overlaps of LFWFs in the 
mixed representation ipn from Eq. fIDSp and flD6p . Note that as the argument of the 
wave functions are ^3^, these equations cannot be expressed as an overlap of the FT wave 
functions $„(crj,6^) in position space. 

However, in the mixed representation, one can write. 



d%^e' 



dxRsi{x, (, 6"' 



dx I dH^e^'^"-'" 



R{xXM 



dxR22ix, C, fe"' 



dx I d%'^e'^'^ ''^ 



p{x, b' 



(D8) 
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dxR22{x, C, b 



C 



dx I d%^e'^ ' R{xX,b' 



dx I d%'^e''^'^-''^p{x,b'^). (D9) 



'0 J JO 

R{x, C,b-^) and R{x, C,b-^) can be obtained in terms of off-diagonal overlaps of LFWFs 
ipni^i, bf) and can be obtained from the above equations. p(x, b'-^) and p(x, b'^) are Soper's 
distributions in the frame C = 0. It can be shown that 

d (foi + ib'^^ 



p{x,b 



-2iM- 



-p{x,b 



(DIO) 



db b 

Eqs. (ID8I) and (1D9I) show the relation between the generalized correlation functions R and 
R with the Soper distribution due to covariance of the form factor. However, the functions 
R and R do not have a probability interpretation, unlike Soper's distribution. 

In the imaginary part of the DVCS amplitude, we have the GPDs integrated with a delta 
function. 



Im[M+^ 



N I dx 
'o 



s{x-C) 



— ^(2^2) 



1-1 



4(l-i)v/W 



E^2^2){xX,t)S{(-x) 



(Dll) 



In terms of the correlation functions defined above, this can be written as 



Im[M^ 



N I d%^e'^ ■''^ I dx 



1-C 

+9{x-C)R22{xX,b^) S{x-0~ 



9i(-x)Rsi{x,C,b^) 



e(C-x)i?3i(x,C,6^) 



+^(x-C)i?22(x,C,n 5{x-0 



dx 







Ci?(x,C,&^)5(x-C) 



is the normalization constant. Integrating over x we obtain 



(D12) 



i-CR{x = cy 

R{x = cy 



(D13) 



Thus, the FT of the imaginary part of the DVCS amplitude with respect to gives both 
R{x, C) and C, &±) where x of the struck parton is now fixed at a; = These, in turn. 
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are related to Soper's distributions p{x,b^) through Eq. (IDSP and (lD9p . This is a mixed 
coordinate and momentum space representation. The above relation can be generalized to a 
hadron in a model independent way. Introducing the complete 3D spatial amplitude p{a, b^) 
at fixed light front time r, we can write, 

/m[M++] = N j dae-'"^ J e-'^^ '^"^ pi(J,b^). (D14) 

Here a is conjugate to (. Note that as we are at fixed r rather than at fixed time, there is 
no conceptual problem due to Lorentz boosts. 

The physics of the real part of the DVCS amplitude is more involved. However, it is 
related to the imaginary part by a dispersion relation in x. The real part can be expressed 
in terms of the densities p and p as well, however, it contains a principal value (PV) integral 
over X. We can call the result of the PV integral 7(C, Again, after taking a FT in ( we 
obtain the amplitude in full 3D coordinate space. 
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